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Abstract: We consider the application of the DRA method to the case of several 
master integrals in a given sector. We establish a connection between the homoge- 
neous part of dimensional recurrence and maximal unitarity cuts of the corresponding 
integrals: a maximally cut master integral appears to be a solution of the homoge- 
neous part of the dimensional recurrence relation. This observation allows us to 
make a necessary step of the DRA method, the construction of the general solution 
of the homogeneous equation, which, in this coupled system of difference 

equations. 
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1 Introduction 

Recently, a method of evaluating Feynman integrals based on the use of dimensional 
recurrence relations [1] and analytic properties of Feynman integrals as functions of 
space-time dimension d (the DRA method) was introduced [2]. It was successfully 
applied in a series of calculations [3-9] where master integrals for various families of 
Feynman integrals were evaluated exactly in d (in terms of nested sums) and also up 
to high order in e = 2 — d/2 (in terms of the conventional multiple zeta values (MZV), 
using PSLQ). This fast advance of the DRA method was partly due to the availability 
of a number of magnificent tools and methods: IBP reduction tools, in particular, 
FIRE [10], the sector decomposition analysis of singularities implemented in FIESTA 
[11], the application of Mellin-Barnes technique [12-17], the PSLQ algorithm [18]. 
The DRA method provides results in the form of converging (uniformly in d) nested 
sums with factorized summands. Such a form allows one to evaluate many terms 
of the e-expansion with very high precision. This feature of the DRA method was 
demonstrated in the evaluation of master integrals for four-loop massless propagators 
which were previously evaluated in [19] up to transcendentality weight seven. Using 
the results of DRA method it was possible to perform an evaluation up to weight 
twelve [9], and it is certainly possible to go further. 

However, up to now, all applications of the DRA method concerned cases with 
only one master integral with a given set of denominators (in a given sector). The 
reason is that the DRA method requires finding the general solution of the homo- 
geneous part of the dimensional recurrence relation. For the case of several master 
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integrals in one sector this problem becomes very nontrivial. The corresponding ho- 
mogeneous equation has a matrix form and is equivalent to one difference equation of 
order higher than one. One may speculate that this problem is, in a sense, artificial 
and the homogeneous system can be decoupled or, at least, reduced to a triangu- 
lar form by a proper choice of the master integrals (i.e., by passing to some linear 
combinations of the integrals with rational coefficients). In this case the high-order 
difference equation for one master integral should have a hypergeometric-term solu- 
tion, which can be checked by the Petkovsek's algorithm Hyper [20]. Unfortunately, 
in real-life examples, the homogeneous equation appears to have no hypergeometric 
or d'Alembertian solutions. So, taken as a separate mathematical problem, finding 
the solution of the homogeneous equation for the case of several master integrals 
cannot be performed in a systematic way. Therefore, when constructing the homo- 
geneous solution, one has to rely on some additional methods. The goal of this paper 
is to present a method to find the homogeneous solution using unitarity cuts. The 
idea is very simple and yet appears to be very useful. 

The paper is organized as follows. In the next Section we introduce our nota- 
tion. In Section 3 we show that the maximal cut of the master integral satisfies the 
homogeneous part of the dimensional recurrence relation and this property gives a 
practical tool of finding a solution of the homogeneous part of the dimensional re- 
currence relations. In Section 4, we illustrate our technique on the evaluation of two 
master integrals (called Iu and J15 in [21]) for the three- loop static quark potential 
[21-24]. We reproduce the results presented in [21] and obtained with the help of 
the Mellin-Barnes representation [12-17] and obtain one more term in e-expansion 
(weight seven). 

2 General setup 

Let us suppose that we are interested in the evaluation of an L-loop Feynman integral 
depending on E linearly independent external momenta pi, ■ ■ ■ ,Pe- There are N = 
L(L + l)/2 + LE scalar products involving the loop momenta If 



where e a = ±1 and v = d/2 is a convenient variable. The quantities e a i0 = ±i0 
determine the infinitesimal shifts of the denominators poles. The scalar functions 
D a are linear polynomials with respect to s^. The functions D a are assumed to 
be linearly independent and to form a complete basis in the sense that any non- 
zero linear combination of them depends on the loop momenta, and any Sik can be 



Sik = h-qk, i = l,...,L, k = l,...,L + E, 
where qi,...,L = h,...,L, <?l+i,...,l+e = Pi,...,e- The integral has the form 
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expressed in terms of D a . The indices n a are assumed to be integer, and if n a > 
we say that the integral has a denominator D a . The integrals having the same set 
of denominators form a sector. 

The dimension shifting relation can be written in two equivalent forms [4, 25]: 

J {y - 1; n) = Q (A u . . . , A N ) J {u; n) , (2.3) 

or 

J (u + 1; n) = P (B u . . . , B N ) J (u; n) , (2.4) 

where Q (A±, . . . ,An) and P (Bx, . . . , Bn) are some polynomials. The operators A a 
and B a act as follows: 

A a J {y\ m, . . . , n a , . . . n N ) = n a J (v; m, . . . , n a + 1, . . . n N ) , 

fi Q J (i/; m, . . . , n a , . . . n N ) = J (v; m, . . . , n a - 1, . . . n N ) . (2.5) 

In order to obtain the dimensional recurrence relation for some master integral 
Ji (u) = J {v; ni), we have to plug it in Eq. (2.3) and reduce the right-hand side using 
IBP identities. Observe that the integrals appearing on the right-hand side of Eq. 
(2.3) belong to the same sector as J (is; n) or simpler (lower) sectors. Therefore, the 
result of the IBP reduction is also a linear combination of master integrals belonging 
to the same, or simpler, sectors. Therefore, if there are no other master integrals in 
the same sector as J±, the general form of the dimensional recurrence relation is 

J 1 {u+l)=C(u)J 1 (u) + R(u) ) (2.6) 

where R (y) contains only simpler master integrals, and C (u) is a rational function. 
Naturally, the dimensional recurrence relations for simpler master integrals do not 
depend on J\ . The homogeneous part of this equation can be easily solved in terms 
of T-functions. The situation is different if there is more than one master integral 
in a given sector. In this case we will refer to the column of master integrals in a 
given sector as a multicomponent master integral (MMI). The dimensional recurrence 
relations for MMI form a coupled system of equations which can be written in the 
matrix notation as 

J(j/ + 1) = C(v) J(i/) + R», (2.7) 

(Ji\ 

where J = is an MMI and C (u) is a matrix with rational elements. 

In order to apply the DRA method, we have to find a general solution of the 
homogeneous equation (v + 1) = C (y) Jh {v). This system of difference equations 
can be reduced to one difference equation of fc-th order, for example, for Ji^. In 
particular, for k = 2, we have 

Ji,h (y + 2) + C 1 {u) J 1)h {v + 1) + C 2 {y) J x , h (f ) = 0, (2.8) 
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where C-\_ i(y) are rational functions expressed via matrix elements of C(u). In 
general, difference equations of a high order cannot be solved analytically. There is, 
however, a possibility to check whether the equation has a solution in the form of 
a hypergeometric term (i.e., such a solution / [v) that / (u + 1) // (u) is a rational 
function). This possibility is based on the Petkovsek's algorithm Hyper [20]. In fact, 
the (non-) existence of a hypergeometric-term solution allows one to claim also the 
(non-) existence of a more general solution — d'Alembertian solution. Unfortunately, 
the application of the Hyper algorithm to real-life examples (in particular, to the one 
considered in Section 4) proves that such solutions do not exist. Therefore, solving 
homogeneous matrix difference equations is a very nontrivial mathematical problem. 

3 Cut integrals 

Let us now improve our notation indicating in the arguments of J also the signs of 
infinitesimal imaginary shifts in Eq. (2.2): J [v\ n) — > J (u; n; e). Below we omit the 
first argument v where it does not lead to confusion. 
Then, an a- cut integral can be defined as 

A a J (n; e) = J (n; . . . , e a , . . .) - J (n; . . . , -e Q , . . .) . (3.1) 

Similarly, we can define an integral cut over several lines: A{ aij ... an } J (n; e) = A ai . . . A Qn J (n; e). 
For non-positive n a the infinitesimal shifts do not change the integral, so, 

A Q J(n; e) = if n a < 0. 

In contrast, for the positive n a the cut integral A Q J(n;e) is not zero and can be 
obtained from J (n; e) by the replacement 

(— l) na 

[D a + e a i0}- na [D a + e a i0]- na - [D a - e a i0]- n " = 2Txie a \-^5^-V (D a ) , 

r (n a ) 

(3.2) 

where 5^ (x) = -^5 (x) denotes the n-th derivative of Dirac's 5 function. For n a = 1 
this prescription reduces to the well-known replacement (D a + iO)^ 1 — > —2iri5 (D a ). 

It is clear that the IBP identities are not sensitive to the specific choice of e a 
in the sense that the coefficients in these identities do not depend on e. We note, 
however, that the symmetry relations are sensitive to this choice since a symmetry 
that replaces D a — > Dp also replaces e a —> e@. Therefore, temporarily we consider 
the master integrals, identical due to the symmetries, to be different. Then the IBP 
reduction of an integral J (n) is also insensitive to the choice of e a , i.e. 

J(n,e) = J2C t (n)J 4 (e), 

i 

where only master integrals Ji (e) depend on the choice of e, but not the coefficients 
C % (n), which are rational functions of n, d, and external invariants. The master 
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integrals entering the right-hand side either belong to the same sector as J (n, e) or 
to simpler sectors. Applying the A Q operator to this equation we nullify all integrals 
without Da-denominator. Thus, cutting all the denominators of J (n, e) keeps on 
the right-hand side only master integrals of the same sector as J (n, e). 

The cut integrals are also the basic tool of the powerful generalized unitarity 
technique [26, 27] which provides the possibility to construct scattering amplitudes. 
In fact, this strategy of writing an Ansatz as a linear combination of some basic scalar 
integrals and constructing the corresponding coefficient functions is very similar to 
the strategy of solving IBP relations, especially within Baikov's method [28, 29]. 

The dimensional recurrence relations are also insensitive to the choice of e. This 
is obvious already from the fact that, at the derivation of this relation, we never 
needed to specify explicitly the shifts ±z'0 in the denominators. Therefore, restoring 
e-dependence in Eq. (2.7), we obtain 

J {y + 1; e) = C (u) J (v; e) + R (v, e) , (3.3) 

where the inhomogeneous term R (y, e) = C (y) Jj {y; e) includes only integrals 
of the lower sectors. The matrix C (y), as well as C- 7 (y), do no depend on e. Taking 
the cut A{ \ over all the denominators of J, we nullify this term and obtain: 

A { ... } J {y + 1; e) = C (y) A { ... } J (u; e) . (3.4) 

Thus, we arrive at a simple but important observation: the maximal cut A/ x J [y] e) 
of a MMI J [y\ e) is the solution of the homogeneous part of the dimensional recur- 
rence relation for J {y\ e). 

Two remarks are in order. First, (^-functions in a cut integral may be too re- 
strictive to give a non-zero result for a specific choice of the metric signature. This 
becomes obvious in the Euclidean case, where the denominators are always posi- 
tive. But this is also true for Minkovskian metrics as we will see below. Therefore, 
to satisfy all the restrictions imposed by the 5-functions one may have to choose a 
more general metric signature (1,1,...,— 1,...). Second, a cut integral gives only 
one solution of the difference equation, while for a fc-th order equation, there are k 
independent solutions. For k = 2 we can, in principle, find a second solution in an 
algorithmic way by the 'constant variation' method. However, we find it possible, 
and even more convenient, to guess several solutions of the homogeneous equation 
by examining the Mellin-Barnes representation for the cut integral. The guessed so- 
lutions can then be checked to satisfy the homogeneous equation either numerically, 
or strictly, by using Zeilberger's method of creative telescoping [30, 31]. 
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Figure 1. Master integrals F\ and F% 
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Figure 2. Simpler master integrals P51, P53, P% 2 

4 A three-loop example 

Let us evaluate the two master integrals shown in Fig. 1: 

(i7i d ^y 3 d d k d d l d d r 




(-k 2 )(-r 2 )(-(l + q) 2 ) a (-(k - /) 2 )(-(/ - r) 2 )(-v ■ k)(-v ■ r) 



(4-1) 



where a = 1 and 2, the external momentum q is of the form (0, q), v — (1, 0), and 
— iO is implied in all the propagators. 

The simpler master integrals are depicted in Fig. 2. Here we follow the labeling 
of the master integrals applied in our future paper [32]. Moreover, in this labeling, 
F\ = Pji and Fx = P72 but we keep the notation Fi which is more convenient within 
the present paper. The dimensional recurrence relation reads: 



where F (u) 
and C (iA = 



F(u + 1) 



, R» 



C (u) F (u) + R (u 



(4.2) 



R*(v, 



depends on simpler master integrals, 



is a matrix with rational elements. The functions 



C n (u) C 12 (uY 
C21 {y) C 22 (v) 

Cij (y) and R (y) are presented in the Appendix. Observe that although [y] are 
quite cumbersome, the determinant of the matrix C {y) has a simple factorized form: 



det C {y) 



>-2)(4z/-7) 2 (4z/-5)^ 



16(i/ - 1) 5 (2^ - 3) 2 (8z/ - 13)(8f - 11) (81/ - 9) (81/ - 7) 
This seems to be a general situation. 



(4.3) 
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The homogeneous equation reads 



F h (i/ + l) = C(i/)F fc (i/). (4.4) 

The solution of this equation is equivalent to the solution of the second-order differ- 
ence equation for Fi t h 

F 1>h (u + 2) + d (u) F llh + + C 2 (u) F 1>h (u) = 0, (4.5) 

where C\ and C2 are known functions. As we already mentioned earlier, the solution 
seems to be out of reach of the conventional mathematical methods based on the use 
of the Hyper algorithm. 

In order to apply the DRA method, we need to find two fundamental solutions 
of Eq. (4.4), forming a matrix F/, (v) = (F\ (y) , F^ (v)). Then, using the method 
described in Section 5 of Ref. [2], we can find the summing factor S(f), satisfying 
the equation 

§ {y) = § (y + 1) C (U) . (4.6) 

As we explained in the previous section, the maximally cut MMI AF satisfies the 
homogeneous equation (4.4), or, equivalently, AF\ (y) satisfies Eq. (4.5). Observe 
that contracting the lower line of F\ in Fig. 1 we obtain a scaleless integral which is 
zero. Therefore, there is no need to cut this line as this cut nullifies no simpler master 
integrals. In what follows, we also omit the factors —lixi from each cut. Thus, we 
consider Fx and perform the replacements l/{k 2 + iO) — > 5(k 2 ) and l/(v ■ p + iO) — > 
S(v ■ p) = 5(po) for all the propagators apart from 1/ (— (/ + g) 2 ). 

Let us, first, integrate over the loop momenta of the two identical one-loop 
subdiagrams consisting of one static and two usual propagators 

J (0 = / ^W(k 2 W 2 - 2/ ■ fc) , (4.7) 

where S(k ) comes from l/(v ■ k + iO) = l/(k + iO). Here is a subtle point because 
in Minkowskian metrics we might conclude that this integral is zero due to the 
kinematical restrictions. Indeed, in Minkowskian space the first two 5-functions 
result in k — 0, which is incompatible with the last 5-function. Let us instead use 
the metric signature (1, 1, —1, —1, . . .), so that k 2 = k^+kl—k 2 — . . .—k\ = kl+k 2 —k 2 . 
Then a straightforward integration gives 

2 Mid -2) (-l 2 ) d - 4 . . 



7T 



d/2 Q2Vd-3)/2 



where l 2 = — 1\ + P, and Cl(d) = 2n d ^ 2 /T(d/2) is the volume of the unit hypersphere 
in Euclidean d- dimensional space. 
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To take the final integral 

. , > 1 f d d l J(l) 2 . , 

we turn to Euclidean space and separate the two terms in the denominator of 1/ il\ + 
(1 + q) 2 ) introducing a onefold MB representation. The factor corresponds to six 
'time-like' integration variables, two per each loop momenta. 

Then the internal integration is taken straightforwardly and we arrive at the 
following result: 

2 4 ~ 4l T (6 - 3u) 

1 [u) ~ i>-i) 2 r(8-4^-f) 

dz T(-z)r (z + l)T (3u - ^ - z)T (z - 4u + 8)T (z + u - 1) 
2tH T(z + 5-2u) ' 

(4.10) 

It is easy to convert this representation to a linear combination of 3-F 2 hypergeometric 
functions. 

As we mentioned earlier, this gives us only one solution, while a second-order 
equation should have two linearly independent solutions. In order to find both solu- 
tions, let us observe that there are two series of poles from the right of the integration 
contour and three series of poles from the left: 

11 

Zi — n, z 2 = 3v — — + n, 

Z3 = — — — n, z^ = Au — 8 — n, z§ = 1 — v — n, 

where n = 0, 1, ... It turns out that the contribution of any of these series constitute a 
solution of Eq. (4.5). This can be checked either numerically, or using the Zeilberger's 
method of creative telescoping [30, 31]. We assume, of course, that the corresponding 
sums are defined in some region of v where they converge, and then analytically 
continued to the whole v complex plane. As two independent solutions we choose 
the contribution of the series of residues at z\ and Z4. The solutions have the form 



V ^2 4 - 4 T(6 - 3u)T (3i/ - f ) / 8 - 4v ± v - 1 



T(5 - 2z/)I> - l)r (4u - f ) ^ 2 V 5 - 2u, f-3u 
32r(6 - 3u)T(5u - 9)r (| - v) / 8 -4r/,|-z/,4-2z/ 



F iA v ) 2^(8u- 15)r(z/- l) 2 r(2z/-3) sF2 V 10-5*/,f-4z/ 



(4.11) 
(4.12) 



Analytical properties of F± h (u) and F 2 h (y) can be found from the above repre- 
sentation. Conventional series representation of the hypergeometric functions 3F2 in 
Eqs. (4. 11), (4. 12) converges at dtu < 5/2. In order to determine the analytical prop- 
erties of Fl h {y) and F± h (u) in the region dtu ^ 5/2, one has to use the recurrence 
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relation (4.5). It would be more convenient to use the representation in terms of 
series converging uniformly in d. Luckily, both 3F2 in Eqs. (4.11) and (4.12) appear 
to be nearly-poised, and it is possible to transform them to Saalschutzian 4F3, whose 
series converge uniformly in d. Explicit expressions of F^ h (v) and F 2 h (z/) in terms 
of Saalschutzian 4 F 3 are presented in the Appendix. Therefore, the fundamental 

(F 1 (u) F 2 (u)\ 
\ h \ h , where F% h (u) and 
F 2,h i v ) F 2,h \ v )) 
F% h (u) are obtained form the first equation of the system (4.4): 



F* (1/ + 1) - C„ (v) F* (v) ^ 2 , F* {u + 1) - C u (u) F{ 



F ^ (v) = c^7) ' F ™ (v) c a (v) 

(4.13) 

Now, following the recipe formulated in Section 5 of Ref.[2], we obtain the summing 
factor 



S(,)=W MS („^- p r w{) , (4,4, 

where S (u) = 2 2v {v - 2)T(2u - 3) 2 T (4z/ - f ) / (t {2v - T(2 - vf sin(7ri/)) is a 
solution of the equation S (u) = S [y — 1) det C [y) and W iy) is an arbitrary periodic 
matrix. Using Eqs. (4.2) and (4.14) , we obtain the relation 

(SF) {v - 1) = (SF) (v) + S (v - 1) R (v) , (4.15) 

which implies 

(SF)(i/) = W(v) + £ S(v-1)R», (4.16) 

+00 

where W {y) is an arbitrary periodic column-vector and the notation E±oo / iy) in- 
troduced in Ref. [25] means 



s fiy) = -J2f(" + n), 

hoo z — ' 

n=0 

00 

S /(") =X)/(f-n). (4.17) 



n=l 



Now we need to determine W (i/) from the analytical properties of (SF) (1/) 
which depend on our choice of W(z/). In particular, if we choose W (u) = 1, the 
function (SF) has singularities at 1/ = 2, 2±, 2±, 2|, 2|, 2|, 2|, 2|, 2§, 2§ on the stripe 
3?^ G [2, 3). In order to cancel these singularities, we can choose W {v) to be properly 
degenerate (and sometimes completely vanishing) matrix at the points of singulari- 
ties, but we should also try to not spoil the behavior of (SF) at v — > ±200. There- 
fore, it is very useful to eliminate also the explicit and hidden zeros of S, which, at 
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1, are located at the points v = 2|, 2|, 2|, 2|, ±ioo. We finally choose 



1 + c)(l + 2c) / 2 5 (1 - c) (1 - 2c - 4c 2 ) 2 5 1 i±£ T (l - 2c) 



"I -^(l-2c-4c 2 ) ^y(l-2c-4c 2 ) 



(4.18) 

where c = cos {2irv). With this choice of the summing factor, (SF) is holomorphic in 
the stripe dtu G [2, 3) and grows at v — > ±ioo slower than exp (27r \ u\). Taking into 
account the singularities of S+oo S (z/ — 1) R (v), we obtain 

W 0) = , ^ - (it - 2arctan (4^) cos 2 (tt^)) . (4.19) 

Multiplying Eq. (4.16) by § _1 (u), we obtain 

F (u) = (i/) W (v) + S" 1 (z/) E § {v - 1) R (u) . (4.20) 

With quantities § (y) , W (u) , R(z/) determined by Eqs. (4.14), (4.18), (4.19), and 
(A. 2), the above representation (4.20) gives the final result of the DRA method for 

F 2 (v) 

Let us make two remarks about the two terms in this representation of ( z/ )- 
The second term, in fact, does not depend on the explicit form of the summing factor 
§ (v) because 



the MMI F (v) 



is always a finite product of rational matrices. This product can be evaluated recur- 
sively, so that one can organize a numerical evaluation without nested loops. The 
first term can explicitly be written as a combination of fundamental solutions F^ 



and F 2 : 



- 1 (u)W(u) 



Fi, h {y) 
F2,h (y) 



2 5 vr 5 / 2 (2c- 1) (tt - (c + 1) arctan (4^)) 
Fl,h = (2c+l)(l-c)c 2 

4c3 2 ; 2 2 _ C 1 +1 ^-(^ 2 + 2c-l)^ 

F 1)h (i/ + l)-Cn(i/)F 1>h (i/) 



2,h 



C12 {») 



c = cos(2tti/) . (4.21) 
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Now, taking into account that the evaluation of all the nested sums appearing in 
representation (4.20) can be organized in one loop, it is easy to calculate F (z/) with 
high precision and apply the PSLQ algorithm. Then we obtain: 

( 1808C(5) 8C(3) 2 928tt 2 C(3) . /1\ 8 „, , 

+ f -192s 6 + ±±± _ + ^1 + 647r 2 Li4 l\ + _ n 2 ln 4 (2) 

20 4l 2/n s QO 4l , , 428tt 6 1792tt 4 

7r 4 ln 2 (2) - 327T 4 ln(2) 1 

3 v ; v ; 2835 135 

+ (-768Li 4 f - J C(3) - 1287r 2 Li 5 (J) + 512vr 2 Li 4 (-) - 1536s 6 + ln(2) 



2 J \2J 7 

384s 7a 3072s 7fe 4960C(7) 35519vr 2 C(5) 14464C(5) 64C(3) 5 



7 7 21 42 3 3 

32C(3) ln 4 (2) + 372C(5) ln 2 (2) + 32tt 2 C(3) ln 2 (2) 



31457tt 4 C(3) , 7424vr 2 C(3) , , M ,, /QU , 



945 9 

480,^,, fn . 3424vr 6 14336tt 4 16 2l 5 , n , 64 2l 4/n , 

C(3) ln(2) + + — vr 2 ln 5 (2) + — rr 2 ln 4 (2) 

7 SV ; V 1 2835 135 15 V ' 3 V 7 

- -7T 4 ln 3 (2) - — rr 4 ln 2 (2) - — rr 6 ln(2) - 256vr 4 ln(2)) e 2 + 0(e 3 ) , (4.22) 
9 3 315 J 



4 

F 2 (2 - e) *™ -— - 93C(5) - 14vr 2 C(3) - 2rr 4 ln(2) 



-96s 6 + 120C(3) 2 + 327r 2 Li 4 ( £ ) + \v 2 ln 4 (2) - ^tt 4 ln 2 (2) 



3 420 
.384 Lil (I) C (3) - 64*4 (i) + ^ M 2) - ™± - ™* 
32666C(7) 40585tt 2 C(5) , 35047tt 4 C(3) 1fi ., ,_ 4 , , , 2 , 



7 84 630 



16C(3)ln 4 (2) + 186C(5) ln 2 (2) 



+ 167r 2 C(3) ln 2 (2)-^C(3) 2 ln(2)+ Att 2 ln 5 (2)-^r 4 ln 3 (2)-^7r 6 ln(2)) e 2 +0(e 



(4.23) 

where the notation = indicates that the equality holds numerically with at least N 
decimal digits, 



s 6 = C(-5,-l) + C(6), 

s ra = C(-5, 1, 1) + C(-6, 1) + C(-5, 2) + C(-7) , 
s 7b = C(7) + C(5, 2) + C(-6, -1) + C(5, -1, -1) , 
and C( m i) • • • ) m k) are multiple zeta values 

oo ii — 1 ifc-i — 1 / xj. 

«.)=EE''En^' (*■*) 

ii=lia=l «fc=l j'=l *j 
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The terms up to e 1 are in agreement with the previous results [21]. 



5 Conclusion 

We have presented here a method of finding the solution of the homogeneous part of 
dimensional recurrence relations for multicomponent master integrals. The method 
is based on the fact that the maximally cut master integral satisfies this homogeneous 
equation. Strictly speaking, it gives us only one solution, while for a /c-component 
master integral we need k linearly independent ones. However, it appears that in the 
Mellin-Barnes representation of the cut integral each series of poles separately gives 
rise to a solution. For each individual case, this fact can be checked both numerically 
and strictly, using Zeilberger's algorithm of creative telescoping. 

As an application of this technique, we have presented the calculation of the two- 
component master integral ( ^J 1 J given by Eq. (4.20) and entering the three-loop 

V72/ 

static quark potential. Using this result, we have calculated with a high precision 
the e-expansion up to e 2 -terms and applied the PSLQ algorithm to express it in 
terms of conventional constants. Our next natural task is to complete the analytical 
evaluation of all the master integrals for the three-loop static quark potential, i.e. 
to evaluate the last three analytically unknown expansion coefficients entering the 
corresponding result. 
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A Coefficients in the dimensional recurrence relation 

Here we present for completeness the quantities CV,- (v) and Ri (u) entering Eq. (4.2): 
78656i/ 6 -709872z/ 5 + 2652380z/ 4 -5251197z/ 3 + 5809568z/ 2 -3405384i/+826308 



C*21 



4(i/ - 1) 3 (2^ - 3) 2 (8z/ - 13)(8i/ - 11) (8z/ - 9)(8i/ - 7) 
(u - 2) (13120z/ 4 - 70192z/ 3 + 140108z/ 2 - 123689z/ + 40761) 

8(1/ - l) 3 (2z/ - 3) 2 (8z/ - 13)(8z/ - ll)(8z/ - 9)(8z/ - 7) 
3 (13120z/ 6 - 117424z/ 5 + 434716i/ 4 - 851937i/ 3 + 932032z/ 2 - 539672z/+ 129216) 

8(1 - uf(2u - 3) 2 (8z/ - 13)(8z/ - 11) (8f - 9) 
_{u-2) (6592i/ 4 - 34768z/ 3 + 68324z/ 2 - 59303z/ + 19191) 
22 ~ 16(i/ - l) 3 (2z/ - 3) 2 (8z/ - 13)(8z/ - ll){8u - 9) ' ( ' } 
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Ri (k 
R2 iy) 



Pi (u) = -. r -. 128 w -(648214272i/ 10 - 9064230912z/ 9 

1V ; 5(8i/ - 14) 8 (3i/ - 5)3(21/ - 3)2 V 

+ 56911513696i/ 8 - 211292587888z/ 7 + 513701269195z/ 6 - 854608449763z/ 5 

+ 985285600699i/ 4 - 777347268382z/ 3 + 401666882358z/ 2 - 122748402735// 

+ 16847478900) P 51 (u) 

- 16(2z/-3) -(H33568Z/ 6 - 8922240z/ 5 + 29193664i/ 4 

(8u - 13)7(3i/ - 5)3(1/ - 1) V 

- 50834923// 3 + 49690736z/ 2 - 25855817i/ + 5595660)P 53 {v) 
17440z/ 4 - 94208z/ 3 + 189730z/ 2 - 168891^ + 56091 



My - l) 3 (4z/ - 7) (81/ - 11) (81/ - 9) (81/ - 7) 



P 62 M 



192 

P 2 H = 7 -(21590784z/ 10 - 300317184i/ 9 + 1875217824z/ 8 

2V ; (81/ - 14) 7 (3i/ - 5)3(21/ - 3) 2 V 

- 6922120208i/ 7 + 16728915563i/ 6 - 27658192027i/ 5 + 31682083339z/ 4 

- 24828801753z/ 3 + 12740565282// 2 - 3865560708i/ + 526619520)P 51 {y) 

+ - '— -(47232i/ 6 - 368208i/ 5 + 1192698i/ 4 - 2055044z/ 3 

(8z/ - 13) 6 (3i/ - 5)3(1/ - 1) 



+ 1986779i/ 2 - 1021999// + 218560) P 53 (1/) 

4320z/ 4 - 23000z/ 3 + 45592z/ 2 - 39893z/ + 13008 
+ My - l) 3 (4z/ - 7) (81/ - ll)(8i/ - 9) 

The simpler master integrals are 



P 62 (z/) . (A.2) 



7T 2 CSc(7Tt/) CSc(37TZ/)r(z/ - l) 2 

Psi (z/) " f(5~J=5) ' ( } 

7T 3 CSC 2 (7TZ/) CSc(37TI/)r(z/ — l) 3 

53 ^ = r(4-2i/)r(2i/-2) 2 r(4z/-5) ' 

_ 7r 5 / 2 2 8 - 6v (2 cos(2ttz/) - l)r (| - v) r (2z/ - §) T (2i/ - f ) T (f - 4z/) 
-^62 M = 77; „ NT , 77 11 \ . a7 — x 77 — \ 1- 



(2i/ - 3)r (4i/ - f) sin 3 (7rz/) cos(4ttz/) 2 1 + & '(2z/ - 3)r {y - §) 

(15536z/ 5 -75492i/ 4 + 144596z/ 3 -136177z/ 2 + 62875// - 11340) T [y - \) 
XJ L 5 ■ 4-3"(5z/ - 4) (5i/ - 3)r (2u - \) V (f - 4z/) Psi (Z/) ' 

where E-00 is defined in Eq. (4.17). The result for P 62 presented here is found using 
DRA method. 
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B Homogeneus solutions via Saalschutzian 4F3 

Using Eq. (2.4.2.3) from [33], we obtain 



Kk 0") 



tt/2 r(2^-,) 4 F 3 |^ ^_ i/j3 _ i/ ^_ t/ 
cos(3vrz/)(i/ - 2) 2 r (f - 2i/) r (1/ - |) r(2i/ - 4)r (Av - f ) 



vrr(2 - z/)r (21/ - I) 4F3 



z/, 1 - v, v 



v 



3 5 

4' 4' 



2z/-3 



4cos(37rz/)r (I - 2v) Y{y - \) 2 Y (Au - f ) 



Kh iy) 



(B.l) 



25 

■1 



sin(5vrz/)r (f - 6z/) T{2v - 3)r(3 - 



2 2!y - 5 7rsin(vrz/)r(2 - u) 4 F ; 



\,2-v,2-v,\-v 
*L-2v,¥-2v,v-\ 



sin(57rz/)(8z/ - 15)r (| - u) T{u - l)T(2v - 3) 



fB.2) 
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